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Abstract. We consider a system of trapped spinless bosons interacting with 
a repulsive potential and subject to rotation. In the limit of rapid rotation 
and small scattering length, we rigorously show that the ground state energy 
converges to that of a simplified model Hamiltonian with contact interaction 
projected onto the Lowest Landau Level. This effective Hamiltonian models 
the bosonic analogue of the Fractional Quantum Hall Effect (FQHE). For a 
fixed number of particles, we also prove convergence of states; in particular, in 
a certain regime we show convergence towards the bosonic Laughlin wavefunc- 
tion. This is the first rigorous justification of the effective FQHE Hamiltonian 
for rapidly rotating Bose gases. We review previous results on this effective 
Hamiltonian and outline open problems. 



1. Introduction 

A fundamental characteristic of trapped Bose gases is their response to rotation 
[TS1[TT]. When the angular velocity Q becomes large, a transition from a condensed 
regime to a highly correlated, uncondensed, phase is expected. The behavior of the 
system then has certain features similar to the Fractional Quantum Hall Effect that 
is observed in superconductors submitted to a magnetic field. For rotating Bose 
gases, this regime has not yet been observed experimentally, the corresponding 
value of ^ being unattainable at present. Nevertheless, there has been a lot of 
interest in the theoretical understanding of this phenomenon in the literature (see, 
e.g., [11] for a recent review). 

In this paper we present a rigorous study of such a Bose system with a generic 
repulsive interaction potential. We show that in a certain limit, the Hamiltonian 
of the system can be replaced by a simplified effective Hamiltonian with a contact 
interaction, projected onto the Lowest Landau Level (LLL). This simplified model 
has been extensively used in the physics literature [15l [13l [TH [32l [30l [36j [29]. We 
also prove that in the same limit, the true ground state of the system converges 
to the ground state of the model Hamiltonian. In a certain parameter regime, our 
analysis provides a rigorous derivation of the bosonic equivalent of the well-known 
Laughlin state [2"D] . 

Let us consider N interacting spinless bosons submitted to a rotation around 
the x 3 axis and a harmonic trapping potential. Denoting x — (x , x 2 , x 3 ) € R 3 , the 
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Hamiltonian of the system in the rotating frame is given by 



N r i - 12 



/ W a (xj-x k ). 

l<j<k<N 

Here L = i x p is the angular momentum, e*3 = (0,0, 1), m is the mass of the 
bosons, and uj± and u>ti are the trap frequencies. The interaction potential W a is 
assumed to be non-negative, i.e., purely repulsive, and to have scattering length a 
(see [24] for a proper definition of the scattering length). It is natural to introduce 
a fixed potential W with scattering length 1 and write W a (x) — a~ 2 W(x/a). For 
convenience we will assume that the angular velocity is nonnegative. The above 
Hamiltonian acts on the space of permutation-symmetric square-integrable TV-body 
wavefunctions. 

The Hamiltonian ([I]) is stable (i.e., bounded from below) only when f2 < lu±. 
The regime of rapid rotation that will be of special interest to us corresponds to 
the case when fl is very close to the maximal possible speed lu±, i.e., 

uj\ — Q 
uj := < 1 . 

Ul±_ 

To simplify certain expressions, we will work with an isotropic harmonic poten- 
tial, luj_ = wn. Our results would hold equally well when ujj_ ^ u>\\ but u)\\/u)x > 
e > 0. Similarly, a non-harmonic confinement potential in the x 3 direction could 
be used. 

It is convenient to chose units such that m = H = u>± = 1. Introducing A(x) — 
(— a; 2 ,^ 1 ^) and completing the square, our Hamiltonian |T|) can be written as 



N 



J - + ue 3 ■ L 3 



Wa(Xj-Xk). 

l<j<k<N 



(2) 



The kinetic energy term of this Hamiltonian is equivalent to that of a charged 
particle in a constant magnetic field B = V x A. The spectrum of (\p — A(x)\ 2 + 
\x 3 \ 2 )/2 is purely discrete, its eigenvalues being 3(j + l/2) for j = 0, 1, . . . . They are 
all infinitely degenerate. In the definition ^ of our Hamiltonian we have subtracted 
the unimportant ground state energy 3/2 of the kinetic term. 

When the speed of rotation O is not too close to oj± and the Bose gas is sufficiently 
dilute, the ground state of ([2]) is known to exhibit Bose-Einstein condensation, 
with condensate wavefunction described by ground states of the Gross-Pitaevskii 
functional [HI [33] 

e* w - (* (fciaim^ + „e- s i (3 ) 

where g = AirNa. In the limit N — > oo with < lu < 1 and g > fixed this was 
proved in [HCI2]. 

The properties of the Gross-Pitaevskii ground state in the rapidly rotating regime 
uj — > were intensely studied in the literature, both from a numerical [9[[T2l[3] and 
an analytical [5[[T] point of view. As the speed of rotation increases, more and more 
vortices appear and the wavefunction acquires a higher angular momentum. The 
location of these vortices is conveniently studied in the Lowest Landau Level (LLL) 
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approximation where one restricts tp to lie in the kernel of \p— A(x)\ 2 + \x 3 \ 2 — 3. 
This LLL approximation is justified [2] when lu <C 1 and gu <C 1, the number of 
vortices being then proportional to 




The Gross-Pitaevskii functional is expected to be an accurate description of 
the ground state of the many-body system @ provided the number of vortices is 
much smaller than the number of particles in the system, i.e., when 

Within the LLL approximation, this was recently shown in [27] to be indeed the 
case. In terms of the filling factor v = N 2 / (2L tot ) [2], ((U) corresponds to v ^> 1. 

When a/(Nu>) is not small, a completely different regime is expected. Evidence 
of strongly correlated states was found in exact diagonalization studies of small 
systems. As the rotation frequency fl is increased from to its upper limit luj_, the 
ground state encounters a series of transitions between certain values of the angular 
momentum. The behavior of the system is similar to the Fractional Quantum Hall 
Effect in fermionic systems, and usually modeled by an effective Hamiltonian with 
contact interaction in the LLL. The rigorous derivation of this effective Hamiltonian 
for Bose gases with generic repulsive two-body interactions is the main purpose of 
this paper. 

2. Main Results 

2.1. Derivation of the Effective Hamiltonian on the LLL. The ground state 
energy in the bosonic sector is given by 



(5) 



where V/i L 2 (M. 3 ,C) denotes the symmetric tensor product and H^ a was defined 
above in $2$. 

We will compare the ground state energy of the above Hamiltonian @ with the 
simplified model consisting in restricting the wavefunction to the N-body Lowest 
Landau Level (LLL) and replacing W a by A-na times a contact interaction potential. 
The LLL is defined as the ground state eigenspace of the kinetic part of the operator 
© at uj = 0. This subspace of \/f L 2 (R 3 , C) reads 

Sj N ■= I* -..,x N ) = F(x\ +ix 2 1 ,...,x 1 N + ix 2 N )e~ ^ £ L 2 (R 3N ) : 

(zi, . . . , zn) i— > F(zi, . . . , zn) is holomorphic and symmetric > (6) 



where we denote x = (a; 1 , x 2 , x 3 ) € K 3 as before. We will use use the notation 
z = x 1 + ix 2 G C and we will sometimes identify it with {x l ,x 2 ) € M 2 . For any 
W G -f) at , we have by definition 



N 

E 
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The ground state energy of the simplified effective model in the LLL is given by 




(7) 



Note that although it makes no sense to use a delta potential in the original Hilbert 
space, functions in the space S^n are all smooth, hence 



(#,£(X1 -£ 2 )#) = j - j \^(x2,X 2 ,X 3 ,...,X N )\ 2 dX2 



■ ■ ■ dx 



N 



makes perfect sense and defines a bounded selfadjoint operator. As we will discuss 
in the next section, for any u) > and a > there exists a ground state ^> £ fix 
for the problem ([7]). 

We emphasize that ([7]) is not obtained by restricting H^ a to the LLL. For small 
scattering length a such a restriction would lead to a similar expression but with 
the wrong prefactor J W a instead of 4ira in front of the 5-interaction (see Remark 
[T] below). In order to obtain the scattering length, it is important to note that 
the LLL restriction is unphysical on length scales much smaller than the effective 
"magnetic length" , which is 1 in our units. If a <C 1 , the scattering process is 
unaffected by the rotation of the system and hence leads to the scattering length 
as an effective coupling constant. 

Our main result is the following. 

Theorem 1 (Validity of effective LLL model). Let W be a nonnegative radial 
function such that J^ >R W(x) dx < oo for some R > 0, with scattering length 1. 
We define W a ■= a~ 2 W( ■ /a) and 



(i) Upper bound. Assume that rj := k K 3 / 2 a < C 1 one has 

E n {ijJ, a) < E^ ll (lj, a) (l - CaK' 3/2 ^j 



1 



-3/4 



1 



min{l, kN 2 } 
for some universal constant C > 0. 
(ii) Lower bound. Let r = min{l, Ka~ 2 / 3 }. Then 

'Na 1 / 3 



N V J\2\>{2/ v ) 3 / i 



(9) 



E N (cj,a) > ^ LL (w,a) [ l-C 



for some universal constant C > 0. 



a^r 



1 

-in 



W 



ISI^ri/Ba- 8 / 9 



(10) 



What Theorem [T] says is that when k stays away from zero (for instance k > 
fixed) and when a is small enough (depending on the particle number N), then one 
can replace the problem of minimizing H^ a with a generic interaction of scattering 
length a by the study of a simplified Hamiltonian acting on the LLL, with a contact 
interaction of strength 4wa. The latter model has some very specific features that 
we will recall in the next section. 
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For fixed K, the leading order correction in our upper bound is of the order 
aiV 1 / 2 as long as W decays at least as |x| -3-4 / 3 at infinity; it is (aN 1 / 2 ) 36 / 4 if W 
decays as | a?| 3 e for < e < 4/3 instead. The error bounds in the lower bound 
are significantly worse. For fixed k the leading error term is No 1 ' 3 . It remains a 
challenging open problem to derive bounds that display a better N dependence. 
These will require a better understanding of the FQHE regime for large N. One 
would expect that, for k fixed, there exist error bounds that are independent of iV". 

The proof of Theorem Q] uses several previous ideas [24] . The upper bound 
requires a two-scale trial function, as suggested first by Dyson in [16] . in order to 
obtain the scattering length from W a . The fact that the ground state in 9)n for 
(0 is not very well known (contrarily to the condensed Gross-Pitaevskii case) is 
an important obstacle, however. As usual, the lower bound is the hardest part and 
consequently our conditions on a are more restrictive. 

The proof of Theorem [T] will be given in Section 13.21 More general upper and 
lower bounds on Ejy(u),a) are stated in (|28 ]) -([29 ]) and ([43]). respectively. 



2.2. Effective Hamiltonian on the LLL and Convergence of States. The 

result of Theorem[T]can be extended to obtain information not only the ground state 
energy but also on the corresponding eigenfunctions. Before we state our result on 
the convergence of ground states in the limit lo — > and a — > with N fixed, we 
recall in this section several important properties of the effective Hamiltonian on 
the LLL. 

It is convenient to introduce the Bargmann space [7 



B N := i, F 



Z N — > C N holomorphic and symmetric 

/••• / \F(zi, 
Jc Jc 



N < OO 



endowed with the scalar product 



(F,G) B := / dz 



Zl ■ ■ ■ 



dz N F(zi,...,z N )G(zi,...,z N )e ^ 



and its associated norm. It can easily be checked that if F G Bn then the function 



* defined by * = Tr-^e - ^! ^ /2 F e Sj N satisfies ||* 



Ll 2 (r3«) 



IFII 



N 

E 



\p ] -A{x, J )\ 2 + \x 3 \ 



312 



we 3 • Lj 




The delta interaction potential is defined on Bm as follows: 



The prefactor has been chosen to ensure that 



,ZN 



and 
I 

(11) 
(12) 



(F,S 12 F) l 



x)\ 2 dx . 



JR 3 
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It can easily be seen that < S12 < (2tt) 3 / 2 , hence 612 is a bounded self-adjoint 
operator on 82- The model Hamiltonian acting on Sat is defined as 



N 






- Ana <5.y 




l<i<j<N 



(13) 



and its ground state energy equals the LLL energy |[7J) introduced in the previous 
section: 



E N LL (uj,a) :=w£v Blf (H» a ). 



We introduce, for convenience, the notation 

N 

Cn •= y^Zjdzj and A N := ^ 

j = l l<i<j<N 

for the total angular momentum and the contact interaction potential in the LLL, 
respectively. Because of rotation invariance of the interaction these two operators 
commute on Bn, i.e., [Cn, An] = 0. Hence the ground state energy E N hL (ui,a) of 
our Hamiltonian H^j a — ujCn + 47raAjv is obtained by looking at the joint spectrum 
of Ajy and Cn- If we denote by An{L) the lowest eigenvalue of the operator An 
in the sector of total angular momentum L, we get 

E N LL (u,a) = inf {ujL + 4vraA A r(L)}. 

Multiplying any common eigenstate of Cn and An by the center of mass Y2jLi z h 
one sees that ct(Ajv)., , n r , C ct(Ajv)., , ,,. Therefore, L 1— > An(L) is 

nonincreasing. 

A sketch of the general form of the joint spectrum of An and Cn is shown in 
Figure [TJ The possible ground states for H^ a are those whose values of Cn and 
Ajv lie on the so-called yrast curv^ [31) which is the graph of the convex hull of 
L 1 ► An{L). We can write 



\AttkN 2 N 

where, as before k — a/(Nui). Thus the ground state which will be picked by the 
system only depends on the value of k. It jumps from one state to another when 
K is varied. The FQHE regime corresponds to k ~ 1 in which case Cn ~ N 2 and 
Aat ~ N, hence E N LL (ui, a) ~ Na. 

The kernel of A^ is obviously given by 

ker(Ajv)= lF( Zl ,...,z N ) J] ( Zi — Zj) 2 I F holomorphic and symmetric 

I l<i<j<N 



4n the literature the graph of L t— > Ajv(-L) is sometimes called the yrast curve. We keep this 
name for the convex hull which contains all the possible ground states of Hfj a . 
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A N 
N 




Figure 1. General form of the joint spectrum of A^/N and 
C N /N 2 . The dashed curve is the graph of £ ^ A N (£N 2 )/N, 
whereas the solid one is the yrast curve. Points of the joint spec- 
trum lying on the yrast curve are emphasized by thick dots. This 
figure represents only a sketch, numerical studies for the joint spec- 
trum of An and C n can he found in [JUl [351 IS] ■ 



The function which has the lowest angular momentum among these functions (hence 
lies on the yrast curve) is the (bosonic) Laughlin wavefunction 



-FlL(>i> ■■■,z N ) = k N Yi z j) 2 

\<i<j<N 

where kjsi is a normalization factor. It satisfies 

C N F£ U = N(N-1)FZ U . 

The Laughlin wavefunction is the unique ground state of H^j as soon as n > 
Ki(N) := —l/(4:Trdi(N)), where di(N) is the (unknown) left derivative at £ = 
1 - 1/iV of the convex hull of £ ^ A N (N 2 £)/N. The filling factor [U of the 
Laughlin function is given by 

= N 2 = 1 TV^oo 1 

^ " 2(i^ au ,^< au ) - 2(1 - l/N) * 2 ' 

In general, the convex hull An of the function L i— > Ajv(X) is piecewise lin- 
ear and we may define similarly K\(N) > K2(N) > • • • > Kk(N) by the formula 
Kj(N) := -l/(Anrdj(N)) where d\(N) > d 2 (N) > ■ ■ ■ > d k (N) are the successive 
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L 2 {N) 




Figure 2. Value of the angular momentum of the ground state(s) 
of H^ a , depending on the value of n. The Laughlin state is the 
unique ground state for all k > Ki(N). The constant function is 
the unique ground state for all < k < K Jmax (N) . 



left derivatives of the function I i— ► An(N 2 £)/N. To any Kj(N) we can asso- 
ciate a unique total angular momentum Lj(N) which is the highest among states 
lying on the yrast curve and having a left derivative equal to dj(N). The corre- 
sponding eigenspace is easily seen to be the ground state eigenspace of H^ a when 
k e (nj + i(N), Kj(N)). When k = Kj(N), the ground state eigenspace is the one 
containing all states lying on the yrast line with slope dj(N). It does not have a 
unique angular momentum. These statements are illustrated in Figure^ 

The only state having L = is the condensed state F(zi, . . . , zjv) = 1, hence 
Ajv(O) = (27r)~ 3 / 2 iV(iV - l)/2. Also Ajv(1) = (27r)- 3 / 2 iV(iV - l)/2, with unique 
state F(zi, . . . , z^) = J^Li z i- Moreover, it is well known [HI ESI [XS] that 

&n{L) = —±—^ N(N-1-\L) iox2<L<N. (14) 

For the proof, one notes that 812 has eigenvalues and (27r) _3//2 and commutes 
with the relative angular momentum L\i = (z\ — z%) (d Zl — d Z2 ) /2; in fact, 5\2 is 
nonzero only on the subspace where L12 = 0. On symmetric functions of z\ and 
Z2, the smallest non-zero eigenvalue of L12 is 2, hence (27r) 3 / 2 <5i2 > 1 — £12/2. 
Summing over all pairs we get 

The very last term is non-negative, which yields (|14[) as a lower bound. Finally, 
one checks that for 2 < L < N the lower bound is, in fact, an equality for the 
states S(z\ — zqm) • ■ ■ (zl — -^cm), where S denotes symmetrization and zqm ■— 
V 1 V; v , 
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No exact formula for Ajv(L) is known if L > N. For large N and L -C ^V 2 
the yrast line was studied in [27], where it is proved that in this limit the Gross- 
Pitaevskii energy becomes exact. The result in |27] implies that in this regime the 
convex hull of A^(L) is proportional to N 3 /L. 

Very little is known about the yrast curve for L ~ jV 2 , in particular concerning 
lower bounds. Upper bounds have been derived using certain trial states (Pfaffian, 
composite fermions |13[ I35| I34j ) which have been shown numerically to have a 
large overlap with (some of) the true eigenstates of the yrast curve, at least for 
small N. A rigorous understanding of the properties of the true eigenstates is still 
missing, however. It particular, it remains an open problem to investigate whether 
liminfAT^oo di(N) > 0. This would imply that the yrast curve has a discontinuous 
derivative at the Laughlin state. It would also show a certain robustness of the 
Laughlin state, in the sense that this state is the ground state for fixed k > n\ :— 
lim sup tv^oo K i(A), independently of the particle number N. 

This concludes our review of the properties of the effective Hamiltonian (fT3f . To 
state our last result, we will denote by P/v(k) the (finite dimensional) orthogonal 
projector in L 2 (M. 3N ) on the ground eigenspace of the operator (^itk)~ 1 Cn /N 2 + 
An/N, multiplied by ■n~ N l^e~ ^'=i l 3 ^' I 2 . This orthogonal projector is constant 
for all k G (ftj+i {N), Kj(N)). For k > Ki(N), it is just the projector on the iV-body 
Laughlin function 

^)H*Lau)(*Laul 

where 

l<i<j<N 

The following theorem is a rather straightforward consequence of our proof of The- 
orem m 

Theorem 2 (Convergence of States and Fractional Quantum Hall Effect). 

Let k > and N > 2 fixed, and denote by a an y chosen sequence of ground states 
of the Hamiltonian H^ a in L 2 (M. 3 ). Then we have 

lim K„-iV(«X„||=0. (15) 

a— »0 mi i ii 

a/(Nuj)^K 

In particular, if n > Ki(N), one has — * ^La U when a — * and a/(Nu>) — > n, 
up to a correct choice of a phase for . 

Theorem [2] shows that the ground state of any system of N spinless trapped 
bosons with repulsive interactions of scattering length a is, for small enough a and 
rotation speed close to the critical one, well approximated by the ground state of the 
effective Hamiltonian (fT"3|) on the LLL. In particular, the states are highly correlated 
and exhibit a bosonic analogue of the FQHE, with transitions between certain values 
of the angular momentum coinciding with discontinuities of the derivative of the 
yrast curve. 

The proof of Theorem [2] will be given Section 13.31 As will be obvious from the 
method of proof, both Theorems [1] and [5] can be extended to low- lying excited en- 
ergy eigenvalues and their corresponding eigenfunctions as well. The corresponding 
analysis is similar to previous studies of the effective one-dimensional behavior of 
Bose gases in highly elongated traps [39], and we shall not give the details here. 
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3. Proofs 

3.1. Preliminaries. In this subsection we shall gather some useful preliminary 
results which will be needed in the rest of the proof. We recall that Sij was defined 
in (fT"2j) . Similarly the three-body delta interaction can be defined on B 3 as 



, x pw ^ 1 rpfzi+Z 2 +Z 3 Z 1 +Z 2 + Z 3 Z 1 +Z 2 + Z 3 \ 

(5 123 F)(z u z 2 ,z 3 ):=j-^F^ , - , ) ( l<» 

where the prefactor was chosen such that 

-3 



(F,6 l23 F) B = ( [ e' t2 dt) [ dx [ dz\F(z,z,z)\ 2 e 

\JR J JR JC 

It defines a bounded self-adjoint operator on B 3 . In fact, 



2 p -3|z| 2 -3cc 2 



By definition all functions hi $)n are smooth. A way to quantify their regularity 
was provided by Carlen in [10] . We state it in the following lemma. 

Lemma 1 (An inequality of Carlen jTU]). For any p S N, there exists a constant 
C p such that for any holomorphic function f € B\ 



Vz G 



dpf 2 



e-N 2 <C p (l + \z\ 2 n\\ff Bl - (18) 



Proof. The proof is a consequence of the Cauchy-Schwarz inequality and the fol- 
lowing well-known coherent state representation pj 1171 [5] for functions / in the 
Bargmann space B\ : 



f{z)=^~ 1 / # (19) 

Jc 

where (p^(z) — e^ z . □ 



Using the smoothness of functions in the LLL, one can control any interaction 
potential by a contact interaction, up to an error. This was done first in [3]. 

Lemma 2 (Controlling interaction potentials in the LLL). Let F S Bn 

and ^(x u ...,x N ) := tt- n ^F(z 1 , . . . , z N )e~ s "=i l ^ |2/2 € Sj N . Let g e L 1 {M 3 ) be 
non-negative and radial. Then we have 
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< 



l<i^j^k<N 
2 / 



L 2 (R3«) 



l<i^#fc<AT 



F \ + N 3 C 

B N 



(21) 



for a universal constant C > 0. 



Proof. An inequality similar to ([20]) was derived before in [1J. We shall only write 
the proof of (f20l) for iV = 2 (the general case is then obtained by summing over 
pairs). We shall omit the proof of (f2"Tj) which is analogous. 

Defining G[u, v) — F ^) , we have 



(*,5(fi -x 2 )#) 



ff (f 1 -f 2 )|F(z 1 ,z 2 )| 2 e- |Sl1 dnda; 2 
g(\/2|«|) |G(u,u)| 2 e- |u|2 - |, ' |2 dwd«, (22) 



where we have introduced 



1 



S(M) := - / / 9(z, x\ - xl)e-^\ ~^ dx\ dx 



which is obviously radial, i.e., depends only on \z\. 

We split the v integral in (|2"2"|) into two parts, corresponding to |?;| < 1 and 
M > 1, respectively. Consider first the case < 1. Using the fact that G(u, •) is 
even because F is symmetric, a Taylor expansion yields 

G(u, v) = G(u, 0) + v 2 f (1 - t)^(u, tv) dt . 



dv 2 

By the radiality of g the cross term vanishes when integrating over angles, and 
hence 



g(V2\v\) \G(u,v)\ 2 e-M'dv 



M<i 



g(V2\v\) \G(u,0)\ 2 e- M2 dv 



l«l<i 



~g(V2\v\)\v\ 



Q2Q 

(l-t) — (u,tv)dt 



"M'tfo. (23) 



We integrate this identity against e \ u \~ du. The first term becomes 

( / g(z,x 3 )e~W 2 / 2 dx ) (F,6 12 F) B2 < ( [ g) (F,S 12 F) B ^ . 

\J\z\<l/V2 J \Jr3 / 

With the aid of Carlen's inequality (fT5)) , the second term is bounded above by 



C 



[ g(z,x 3 ) \z\Ux \\F\\ 2 B2 < C> j g(x) j^^dx |F|* a 

J\z\<l/V2 " Jl 3 1 + \ X \ 
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Finally, for |u| > 1 we shall use again (fT8|) . this time for p — 0, to conclude that 
du [ dv~g(V2\v\) \G(u,v)\ 2 e-M 2 -^ 2 < C [ g(z,x 3 )dx \\F\ 2 Ba 

J\v\>l J\z\>l/V2 

<C [ g(x)-^—dx ml . 

This completes the proof. □ 

Remark 1. Although we will not need it, we note that (|23[) also yields a lower 
bound: 



(*,5(2i-*2)*> aa > Uj{x)e-^dx 1 (F,8 12 F) B2 , 

where we have used the same notation as in Lemma [H Combined with (|20|) this 
shows that the restriction of the operator e~ 3 g((xi — X2)/s) to the LLL converges 
to (/m3 9)$i2 as e 0. 

It will be important to have some a priori bounds on the ground state energy 
of the effective Hamiltonian (fT3| . The following is certainly not optimal but it has 
the merit of being simple. 

Lemma 3 (Simple bounds on Ej^ h (uj 7 a)). We have, for k = a/(Nui) and 
N > 2, 

aiVCmin j-^,ivj < E^{u,a) < aN min |i,y^iv| . (24) 

Proof. The upper bound is obtained by taking as trial states the Laughlin function 
and the constant function, respectively. For the lower bound, we note that 



l<i<j<N ^ ' 



> c min {loN/2 , 2naN(N — 1)} > — min x a , aN 2 } 

where c = inf ag 2 (L\ + L 2 + 5\ 2 ) > 0. □ 

Except for the prefactor, the upper bound is expected to be sharp. In other 
words, the lower bound should hold without the factor 1/./V multiplying k , for 
an appropriate constant C. This remains an open problem, however. 



3.2. Proof of Theorem [H 

Step 1: Upper Bound. We start by proving the upper bound, using the varia- 
tional principle. The main difficulty is to get the scattering length in front of the 
interaction potential. As suggested first by Dyson in [TB] , this is done by multiply- 
ing a trial state <& of fj n by a correlated function S accounting for the short scale 
structure of the ground state. Compared to previous similar arguments in [2"5ll37| . 
a new complication comes from the that the trial state $ in Sjn is not a simple 
product function, but itself already a (possibly) highly correlated state of which 
little is known. Fortunately, the information that $ G combined with simple 
bounds on Ej^ L (uj, a) will allow to get the desired upper bound. 
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Let Ffj,uj,a G Bn be a normalized ground state for the LLL Hamiltonian a 
defined in (fiU)) . which is a common eigenvector of £n and Ajy- We consider the 
following trial state: 

^AT.w.a := <SVa &N,io,a (25) 

where 

S Nta (xi,...,X N ) := Yl fa(\Xi-Xj\) 

l<i<j<N 

for some < / Q < 1 which will be defined later, and 

AT /A sr^N \&j 1 2 

$N,w,a,{Xl, ■ ■ ■ ,Xn) ■= 1 F N ,u, a (Zl, . . . ,Z A r)e _2 -= 1 = . 

Note that the norm of &N,uj,a € L 2 (M. 3N ) equals the norm of Fn,uj,o, S Bn- We 
write 

JV 

J=l i<j 

where 

jg- §3 x x| 2 |z 3 | 2 3 f 
^ = 2 + — 2 2 + 3 ^ ' 

We shall also use the notation H^ Q = J2j(hu>)j for short. Using the fact that f a is 
real, we can argue as in [371 Eq- (4.64)] to get the identity 

(*N,u,a, H^ N , u , a ) = E / l ^'^ V ' a|2 |^,a| 2 + ^(Sjf^N^a , H^o*JV,a,,«>- 

Using and that Fn^^ is a normalized eigenvector of £jv, we have 

Hw,0® N,u,a = ^{FN,u,a:^NFN,u,a)s N ^N,u,ai 

hence 

3=1 

We compute 



VfcgjV, a = ^/„(|xfc - £»|) _ ' /(l^m-^nl)- 

i^fc '' Efe Xl ' l<m<n<iV 

{m,n}^{z,fc} 

Using < / < 1 we therefore get 

1 N 1 

^l V ^.«| 2 ^ E /a(l^-^l) 2 +2 E /ad^-^D/ad^-^ 
fe=l l<i<j<N i^j^k 
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Wc finally deduce that 

N,uj,a\\ 2 ($jV>,a, Hu,0®N,ui,a) 

+ U N , U ,N,a, f £ i^a) 2 + W a fi] {\x % - gj\) j $ N , u , N , a \ 
\ \l<i<j<N J I 

+ \l<&N,u>,N,aA Y, f a (\ S i- S j\)fa(\ S k-3j\) I $N,u>,N,a)- 

The next step is to bound the terms on the right hand side of the previous 
inequality. An essential tool is the inequality (|20|) of Lemma [2] which relates, on 
the LLL, the interaction of a smooth potential '^2 i< j g(xi — xj) with that of the 
contact interaction with coefficient J R3 g. As we will see, for a correct choice of f a , 
we will have J[(f' a ) 2 + W a fl) ~ 4™, as desired. Let g a := {f' a ) 2 + W a {f a ) 2 . Using 
Lemma [5] as well as the bound (|T7)) . we obtain 

(^N.^.a , H^ a < £& LL (w, a) l^.al 2 



+ ga - 47TO ll-A^al 2 + \jl^N Q f^j ^ (FjV, w , a , Aw F N ^ a ) Blt 

+ CN 2 [ 9a (g)J?fdx + CN 3 ( j f a (\x\)-^dx) 2 . 
J R3 1 + |a;| 4 \J R 3 1 + \x\ 2 J 

For an upper bound, we can use (Fff jU ,a , Ajv Fn,u,o,)b < £']\r LL ( a -'i a )/(4 7ra ) ; hence 



Let us now choose f a . As in [25], [24] we take, for some b > a to be specified later, 

11 if s > b 

where u a is the solution of the scattering equation 

-u'i(s) + W a (s)u a (s) = 

with it a (0) = and lim^oo u' a (s) = 1. Integrating by parts and using that < 
u a {s) < s and < su' a (s) — u a (s) < a [25l El] we see that 

Ana 



2 



I 9a= I [{f a ? + W a f a ]<- 
JR. 3 JR 3 1 



a/b ■ 

By splitting the integral into a part \x\ < {ab 3 ) 1 / 4 and \x\ > (ab 3 ) 1 / 4 , one checks 
that 

/ ga(\x\)-^rjdx < — ^— \a 2 b 3 + ab 4 f W{x)dx) . 



l + \x\ 4 ~ 1 -a/b 



\x\>(b/a) 3 / 4 
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Note that for a -C b the second term in the last bracket is small compared to the 
first one if W decays faster than | a?| 3 4 / 3 at infinity. We further have 



R3 



,/ < 4nab 
Ja ~ 1 - a 



/b 



and hence 



dx < 



Anab 3 
1 - a/b 



Let us assume, for simplicity, that b > 2a. Then 

, H^ a * N , u , a ) < E^ ll (uj, a ) (l + C [~ + Nab 2 ] ) 

+ CN 2 a 2 b 3 I 1 + Nb 3 + 

\ a J\x\>(b/a) 3 / 4 

for some constant C > 0. To bound the last term relative to the first one, we can 
use the lower bound of Lemma [3] to conclude that 



W{x)dx 



(^N.^.a , H^ a *a^, q ) < E^ ll {uj, a)(l + C 



Nab 2 



nN 2 ab 3 
in{l,K7V 2 } 



1 + Nb 3 + 



a J\S\>(b/a) 3 /' t 



W(x)da 



(28) 



It remains to derive a lower bound on ||'\D r jv',ci>,a|| < Arguing as in :26j we can bound 



II* 



[ ■■■[ n /-a 



x, - xA) 2 \$ 



JV,C 



l<i<j<N 



l<i<j<Af 



(l-/ a 2 )(|x l -f,|)|$ N , w , a | ; 



>l-CN 2 (l-/ a 2 )(|f|) 



1 + \x\ 4 

(1 — /a) ) {FN.u.a, &.NFN,uj,a) , 



where we have used again Lemma [2j For our choice of f a in (f27|) , it is easy to see 
ESI El that 



/ (1 ~ ft) < ^ab 2 
Jr 3 



and hence 



(l-/ 2 )(l^l): 



■ dx < Airab 



l + \x\ 4 

Finally, we can use Lemma[3]to bound (-F/v^.a, &NFN,uj,a)i3 N < Ejf jL (w, a)/(47ra) < 
(2tt)' 3 / 2 Nk- 1 . This yields 



W^NuaW > l-CNab z [- + Nb z 



(29) 



for an appropriate constant C > 0. 

Combining (J5SJ and (EI]), the choice 6 = 2k" 1 / 4 7V" 1 / 2 leads to the desired 
inequality ©. 
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Remark 2. If we had a lower bound 

E^ LL (uj,a) > cNan^ 1 , 

as is expected for k > N^ 1 , the upper bound of Theorem^ could be somewhat 
improved. In %2B\) min{l,K7V 2 } could be replaced by N min{ 1 , k N} in the denom- 
inator, and the optimal choice of b would then be b = (kN)^ 1 ^ 3 . The main error 
term would then be of the order are 1 ' 3 jV 1 ' 3 instead of an 

1/4^1/2 _ 

Step 2: Lower Bound. As a first step, we shall replace W a by the finite range 
potential W a ,R = W a \{\x\ < Rq) for some Ro to be chosen later. Since W a is 
assumed to be non-negative, this is legitimate for a lower bound. We denote by 
a(R ) the scattering length of W a ,R - If -^o is large enough compared to a, we will 
have a(i?o) — a - Indeed, let us recall that [53] 

4™ > 4 7 ra( J R ) = / W aM J aMo > [ W a , R J a >4na- f W a (30) 

JR 3 Jm. 3 J\x\>Ro 

where f a < f a .Ro are the solutions of the zero-scattering equations corresponding 
to W a and W a ^ , respectively. 

We continue with a lemma inspired by a method of Dyson [TO] . The key idea is 
to replace the "hard" interaction potential W a by a softer one using parts of the 
kinetic energy, with this softer potential being close to inaS when projected to the 
LLL. Note that this step is essential, it is not possible to project the original W a 
to the LLL level. This would also lead to a 6 interaction for small a, but with 
the wrong coupling constant J W a instead of 4ira (see Remark [T]). Compared to 
previous studies [551 1231 12H 13H] where a similar strategy has been applied, the 
main new difficulty comes from the fact that the effective kinetic energy ho = 
(\p— e*3 x x\ 2 + |a; 3 | 2 — 3)/2 is not positive locally, i.e., on a domain with Neumann 
boundary conditions, but is positive only on the whole of R 3 . To circumvent this 
problem, we shall rewrite (\l/,/io^') for ^ S L 2 (R 3 ) as 

h V) = \ ( e-l*l 2 (|cW(f)| 2 + \d-M£)\ 2 ) dx (31) 

Z JR 3 

with ip(x) = e' a: ' 2 / 2x I'(a;) and d s — d x i + id x 2. The integrand on the right side 
is now positive but contains no derivatives with respect to z and is hence weaker 
than |V-0| 2 . Nevertheless we shall show in the next lemma that it is still strong 
enough to accomplish the goal of replacing W a by a softer potential for a lower 
bound. The resulting "potential" turns out not be a potential in the usual sense of 
a multiplication operator, but rather is a non-local operator which has the property 
that its projection to the LLL is proportional to a 5-function, however. 

Lemma 4 (Dyson-type inequality). Let y = (s, y 3 ) £ R 3 . For R > Rq, we have 
for all ip 

f e~W 2 (\d x 3*P(x!)\ 2 + \d- z ^{x)\ 2 + W a (x - y)\ij(x)\ 2 ) dx 

J\x-y\<R 

>4™(i? )e-^ 3 l + ^ 2+ l s l 2 



4ttR 2 



e sz %l){x) dx 



\x-y\=R 



(32) 
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Note that if ip € Bi, one has 

4ltR l J\S-y\=R 

and the right side of (l32|) equals 47ra(i?o)e^' J ' |^(2/)| 2 which is precisely 

4:TTa(Ro)(ip,5 y 'ip} 13i when i? = 0. 

Proof of Lemma^ Let <?(x) = e~ sz i(j(x + y). Using that W a > W a ,Ro, we have to 
show that 



\x\<R 



N ) (|^3 S (f)| 2 + |9 2ff (f)| 2 + ^ a , fl0 (f)| 5 (f)| 2 ) 



> 47ra(i? ) e" 



4tT.R 2 



5(3;) <ix 



|af|=fi 



(33) 



Since |z| 2 + (x 3 + y 3 ) 2 < (|?/ 3 | + R) 2 in the integrand on the left, this will follow if 
we can show that 



\x\<R. 



(\d x 3g(x)\ 2 + \d z g(x)\ 2 + W a . Ro {x)\g{x)\ 2 ) dx 



> 4Tra(R ) 



4irR 2 



\S\=R 



g{x) dx 



(34) 



Let now f a ,R be the solution of the zero-energy scattering equation —Af aRo + 
W a ,R fa,R = 0, subject to the normalization Mm^^,^ f a ,R {x) = 1. Since f a ,R Q is 
real-valued, Idg/a^J 2 = \d x if a .Ro? + \d x ^ f a .,R \ 2 , and hence 

{\d x *fa,R {x)\ 2 + \d- z f a ,R, a (S)\ 2 + W a , Ro (x) \ f a>Ro (x ) | 2 ) dx < 4ira(R Q ) . 

\x\<R 

The Cauchy-Schwarz inequality implies that 

( \d x3 g(x)\ 2 + \d z g(x)\ 2 + W a , Ro (x)\g(x)\ 2 ) dx 

1 



\x\<R 



> 



47ra(i? ) 



\S\<R 



(d x 3f a ,R d x 3g + d z f a ,Rod E g + W^R f a ,R a g) 



(35) 



Using partial integration, the zero-energy scattering equation as well as the fact 
that \Vf a ,Ro(x)\ =a/\x\ for \x\ > R this yields (JS}- □ 



As an immediate corollary, we see that for any non-negative function p supported 
on [Rq,R] with p < 1, 



\x-y\<R 



e-W 2 (\d x .^(x)\ 2 + \d-Mx)\ 2 + W a {x ~ y)\^{x)\ 2 ) dx 



> 47ra(i? ) e- (|y / dr p{r) 

JRq 



Airr 2 



\x-y\=r 



(36) 



We shall apply this inequality to the Hamiltonian H^f a , for each particle separately, 
considering the other N— 1 particles as fixed. Consider first particle one, and assume 
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that all particles k > 2 are located at a distance > 2R from each other, i.e., that 
\xk —Xf\ > 2R for all k,£ = 2, . . . , N. Then we get, for all functions F(xi, . . . ,xn), 



d Xl e~^ 2 ^d^Fixu x N )\ 2 + \d- Zl F(x u x N )\ 2 

N 

+ ^2w a (x J -x 1 )\F(x 1 ,... 1 x N \ 



J=2 



>4wa(Ro)^e-^+ R ^+M 2 [ R drp(r) / 

~^ JRo 47rr J\S 



e ZjZl F(xi, . . . ,XN)dx-i 



'\&i—2j\=r 

In general, we can get the same bound if we only retain on the right side the Xj's 
for j = 2, . . . , N which are at a distance > 2R from all the others. Using and 
([3"Tj) . we conclude that, for any < 9 < 1, 

N 

H^ a >Y t {9{h ) j +^-L) i )+^a{Bo){l-B) £ [% , (37) 

j=l l<i<j<N 

where the potential t/y is defined as 

(9,U 12 9) = Jdx a ~'J dx N e-\^ 2 e-^+ R ^+^ 2 e-^^ 2 x 



N 



Y[x\s 

k=3 



k -S 2 \>2R. 



dr p(r) 



Aitr 2 



1 F(x\ , . . .,X N )dX! 



(38) 



with F(fi,...,xjv) = 7r w / 4 *(f ^ 2 )ntie |53 ' |2/2 and f a = («i,as?), af 2 = 

(22,^1)- 

The new potential 5Zi<j ^ s a complicated A^-body term which has the ad- 
vantage of being bounded, however. To be precise, the following bound holds. 



Lemma 5. We have 



| ?7i2|| < sup 



p(r) 



(39) 



r 47rr 2 

Proof. Applying the Cauchy-Schwarz inequality to the x\ integration, we see that 



dxi e Z2Z1 i(,> 



\xi—X2\=r 



< 



\xi — X2\=r 



dxie- |5l|2 |V>| 2 



d Xl e lSll2 - mcZ2Zl 



\xi-x 2 \=r 



It is easy to check that -(|x§| + R) 2 + |z 2 | 2 + |^i| 2 - 2Rez 2 zi < for \x\ — a?a| < R. 
Hence 



*(ari, . . . ,x N )\ 2 ^- dx! ■■■dx N 

47r xi - aca r 



(tf.Efo*) < 

This proves the claim. 

In order to minimize the right side of (|39p . we shall choose 

o 2 

3r" 



p(r) 



R 3 - R% 



for R < r < R. 



□ 



(40) 



Note that J R p(r)dr = 1 
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We shall now apply a standard perturbation theory argument to our Hamil- 
tonian. Let P denote the orthogonal projection onto S)n (the LLL for all the N 
particles), and let Q = 1 — P. Let 

N 

A = ^2(e(h ) j +u;(e 3 -L) 

3=1 

and 

B = 4Tra{Ro)(l - 9)J2Uij 

so that the right side of ([37]) equals A + B. We have A = PAP + QAQ. Since B 
is positive we can use the Cauchy-Schwarz inequality to get the lower bound 

B > (1 - S)PBP + (l - 5- 1 ) QBQ 

for any < 8 < 1 . Using Lemma [5] with the choice (|40| for p we see that 

6iraN(N - 1) 



QBQ < Q ■ 



(Recall that a(R ) < a.) Moreover, for 9 > co, 

QAQ >Q (9-u) . 

We therefore conclude that 

A + B > (1 - S)P(A + B)P + Q(e-w- 
In particular, 

6naN 2 



E N {u, a) > min <^ (1 - S) inf a(A + B) \% N ,0-u 



and it remains to study A + B restricted to at • 

For &(xi, . . . , xn) = -k~ n / 4 YljLi e~\ Xj \ 2 / 2 F(zi, . . . , zjv) a bosonic function h 
f) n , we have 

(*,A*) =oj(F,C n F) Bn 

by (HJ), and 



(*,B*) = 27ra(i?o)iV(iV-l)(l-fl) / dx 2 -(l*al+«) 

X j Wx\x ] -x 2 \>2RdXj\^(x 2l X2 1 X il . . . ,x N )\ 2 (42) 

where we have used that F is analytic in z\ and z 2 . For a lower bound, we use 

~\\ X\x J -x 2 \>2R > 1 — X|x 3 -5 2 |<2i?, • 
J>3 J>3 



Letting 



t R = -Z=l dte-^ 2 >1--^. 
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we finally get 

(%B9) > 47ra(i? )(l - 0)t R ( F, I > % I F 



\ \l<i<i<iV / 



-27ra(i? )A r2 (A r -l)(l-6') / cfcr 2 • ■• J dx N X\x 3 -x 2 \<2r\^(x2,x 2 ,^3, ■ ■ -,x N ) 
By Carlen's inequality (|18[) with p = and 5? 2 , a?4, ■ ■ ■ , xjy fixed, we have 



(43) 



|*(x2,f 2 ,X 3 , . . . ,Xjv)| 2 < C / |^(f 2 ,^2,^3,---,^iv)| 2 ^ 



therefore 



/ ^3X|J 3 -5a|<2Ji|*(^2, ^2,^3, ■■• ,X N )\ 2 < CR 3 / df 3 |* (x 2 , $2, 

Jr. 3 Jm 3 



3 ■ 



x 3 , . . . ,x N )\ 2 



We conclude that 



(a+b) c-«)C,-^).w ^„ i4 



Our final inequality is therefore 

P , w • f (1- <?)(!- 0) (tR-C&N) a(Ro) fU-L , . 
Fn(lu, a) > mm < i% (w, a) , 

I ° 

67raiV 2 



^3-^3) 



(44) 



We now optimize constants. Recall from that _E^ LL (w,a) < k 1 aN. We 
choose i? = a 1/9 , # 3 = 2i?n, = lo + 2a 1 / 3 N and 5 = Qnc^^N such that 



6iraN 2 



8(R 3 - R 3 ^ 



Na 1 ^ 3 . 



'"0 

Assuming that k > a 2, 1 3 this expression is greater than £7j^ LL (w, a). Recalling (|3CT 
the final result is then 



E N (u,a) > E]f h (uj,a) \ l-uj-— / W{x)dx-C 



l\S\>a~»/° 

■ (45) 

If k < a 2 / 3 , the choice R = k 1 / 6 , i? 3 = 2R? h 9 = to + 2NGK- 1 and S = GttNk 1 ^ 
yields the desired bound. This completes the proof of Theorem [TJ □ 

3.3. Proof of Theorem [2j The proof is a simple consequence of the bounds in 
the previous subsection, together with the right choice of parameters. If we choose 
Rq, R and S as above, but 9 to be bigger than the previous choice by an amount 
9' > 0, we conclude that 



H» a >[l-9'-C 



Nal/3 , 1/9 ' 

h a ' r 



1 

47T 



w 1 n*< a n 



N 

vv j ii n 

|£|>r 1 /6 a -8/9 / 

+ (0' + ^ LL (a,,a))Q (46) 
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where IT = n^e^^ 1 l^'l 2 / 2 p denotes the projection from L 2 (M. 3N ) onto B N , and 



r = min{l, na~ 2/3 } as in the statement of Theorem^ Let ITi := n N / 4 e 2 ~^j= 1 /2 P n (k) 
denote the projection from L 2 (M. 3N ) onto the ground eigenspace of the operator 
Cn/N 2 + 4ttkAn/N, and 7at(k) = (A 2 — Ai)/Ai > where Ai and A 2 are the first 
and second eigenvalues of this operator, respectively. Introducing II2 = II — Hi, we 
have 



njfl* Hi + I%M a Il 2 > E^ h (co, a) P n (k) 



Note that 



v ; T1 , / £n . a A N 



n 2 



> N 2 ui 1 



Nujk 



1 



II. 



£ 

N 2 



Attk- 



N 



Nujk 



> ^ LL (w,a)(l + 7jv ( K )) (l 
For a, w, and k — a/(Nu>) small enough, we obtain 

-Na 1 / 3 



Nujk 



H^ a >E^ L (uj,a) 



1 



C 



4tt 



^ LL (ca)HMn;n 2 



Combined with the upper bound ([9]), this clearly yields ||n 2 'I , 1 ^ a |j 
0, as was claimed. 



Oand Q* 



N II 
u.a 



□ 



Acknowledgments. The authors would like to thank Elliott H. Lieb for fruitful 
discussions. Financial support from the ANR project ACCQuaRel of the French 
ministry of research (M.L) and the U.S. National Science Foundation under grant 
No. PHY-0652356 (R.S.) is gratefully acknowledged. 



[l 
[2: 
[3: 

[4: 

[5 
[6. 

[7; 

[9 
[10 



References 

A. Aftalion and X. Blanc, Vortex lattices in rotating Bose-Einstein condensates, SIAM 
Journal on Mathematical Analysis, 38 (2006), pp. 874-893. 

A. Aftalion and X. Blanc, Reduced energy Junctionals for a three-dimensional fast rotating 

Bose-Einstein condensates, Annales Henri Poincare, 25 (2008), pp. 339-355. 

A. Aftalion, X. Blanc, and J. Dalibard, Vortex patterns in a fast rotating Bose-Einstein 

condensate, Phys. Rev. A, 71 (2005), p. 023611. 

A. Aftalion, X. Blanc, and M. Lewin, Unpublished. 

A. Aftalion, X. Blanc, and F. Nier, Lowest Landau level functional and Bargmann spaces 
for Bose-Einstein condensates, J. Funct. Anal., 241 (2006), pp. 661-702. 
L. O. Baksmaty, C. Yannouleas, and U. Landman, Rapidly rotating boson molecules with 
long- or short-range repulsion: An exact diagonalization study, Phys. Rev. A, 75 (2007), 
p. 023620. 

V. Bargmann, On the representations of the rotation group, Rev. Mod. Phys., 34 (1962), 
p. 829. 

G. Bertsch and T. Papenbrock, Yrast line for weakly interacting trapped bosons, Phys. 
Rev. Lett., 83 (1999), pp. 5412-5414. 

D. A. Butts and D. S. Rokhsar., Predicted signatures of rotating Bose-Einstein conden- 
sates, Nature, 397 (1999), pp. 327-329. 

E. Carlen, Some integral identities and inegualities for entire functions and their application 
to the coherent state transform, J. Funct. Anal., 97 (1991), pp. 231-249. 



22 M. LEWIN AND R. SEIRINGER 

[11] N. R. COOPER, Rapidly rotating atomic gases, Advances in Physics, 57 (2008), pp. 539—616. 
[12] N. R. Cooper, S. Komineas, and N. Read, Vortex lattices in the lowest Landau level for 

confined Bose-Einstein condensates, Phys. Rev. A, 70 (2004), p. 033604. 
[13] N. R. Cooper and N. K. Wilkin, Composite fermion description of rotating Bose-Einstein 

condensates, Phys. Rev. B, 60 (1999), pp. R16279-R16282. 
[14] N. R. Cooper, N. K. Wilkin, and J. M. F. Gunn, Quantum phases of vortices in rotating 

Bose-Einstein condensates, Phys. Rev. Lett., 87 (2001), p. 120405. 
[15] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, Theory of Bose-Einstein 

condensation in trapped gases, Rev. Mod. Phys., 71 (1999), pp. 463—512. 
[16] F. J. Dyson, Ground-state energy of a hard-sphere gas, Phys. Rev., 106 (1957), pp. 20-26. 
[17] S. Girvin and T. Jach, Formalism for the quantum Hall effect: Hilbert space of analytic 

functions, Phys. Rev. B, 29 (1984), pp. 5617-5625. 
[18] E. GROSS, Structure of a quantized vortex in boson systems, Nuovo Cimento, 20 (1961), 

pp. 454-477. 

[19] M. Hussein and O. Vorov, Generalized yrast states of a Bose-Einstein condensate in a 

harmonic trap for a universality class of interactions, Phys. Rev. A, 65 (2002), p. 035603. 
[20] R. B. LAUGHLIN, Anomalous quantum hall effect: An incompressible quantum fluid with 

fractionally charged excitations, Phys. Rev. Lett., 50 (1983), pp. 1395-1398. 
[21] E. H. Lieb AND R. Seiringer, Proof of Bose-Einstein Condensation for Dilute Trapped 

Gases, Phys. Rev. Lett., 88 (2002), p. 170409. 
[22] , Derivation of the Gross- Pitaevskii equation for rotating Bose gases, Comm. Math. 

Phys., 264 (2006), pp. 505-537. 
[23] E. H. Lieb, R. Seiringer, and J. P. Solove.i, Ground-state energy of the low-density Fermi 

gas, Phys. Rev. A, 71 (2005), p. 053605. 
[24] E. H. Lieb, R. Seiringer, J. P. Solovej, and J. Yngvason, The mathematics of the Bose 

gas and its condensation, Oberwolfach Seminars, Birkhauser, 2005. 
[25] E. H. Lieb, R. Seiringer, and J. Yngvason, Bosons in a trap: A rigorous derivation of 

the Gross- Pitaevskii energy functional, Phys. Rev. A, 61 (2000), p. 043602. 
[26] , One- dimensional behavior of dilute, trapped Bose gases, Commun. Math. Phys., 244 

(2004), pp. 347-393. 

[27] , The Yrast line of a rapidly rotating Bose gas: the Gross- Pitaevskii regime. Preprint 

arXiv: 0904. 1750, Phys. Rev. A (in press), 2009. 
[28] E. H. Lieb and J. YNGVASON, Ground state energy of the low density Bose gas, Phys. Rev. 

Lett., 80 (1998), pp. 2504-2507. 
[29] S. Mashkevich, S. Matveenko, and S. Ouvry, Exact results for the spectra of bosons and 

fermions with contact interaction, Nuclear Physics B, 763 (2007), pp. 431-444. 
[30] A. G. MORRIS and D. L. Feder., Validity of the Lowest Landau Level approximation for 

rotating Bose gases, Physical Review A, 74 (2006), p. 033605. 
[31] B. MOTTELSON, Yrast spectra of weakly interacting Bose-Einstein condensates, Phys. Rev. 

Lett., 83 (1999), pp. 2695-2698. 
[32] T. Papenbrock and G. F. Bertsch, Rotational spectra of weakly interacting Bose-Einstein 

condensates, Phys. Rev. A, 63 (2001), p. 023616. 
[33] L. P. Pitaevskii, Vortex lines in an imperfect bose gas, Zh. Eksper. Teor. fiz., 40 (1961), 

pp. 646-651. 

[34] N. Regnault, C. C. Chang, T. Jolicoeur, and J. K. Jain, Composite fermion theory of 
rapidly rotating two-dimensional bosons, Journal of Physics B, 39 (2006), pp. S89-S99. 

[35] N. Regnault and T. Jolicoeur, Quantum hall fractions for spinless bosons, Phys. Rev. B, 
69 (2004), p. 235309. 

[36] N. Regnault and T. Jolicoeur, Parafermionic states in rotating Bose-Einstein conden- 
sates, Phys. Rev. B, 76 (2007), p. 235324. 

[37] R. Seiringer, Ground state asymptotics of a dilute, rotating gas, J. Phys. A, 36 (2003), 
pp. 9755-9778. 

[38] R. Seiringer and J. Yin, Ground state energy of the low density Hubbard model, J. Stat. 
Phys., 131 (2008), pp. 1139-1154. 

[39] , The Lieb-Liniger model as a limit of dilute bosons in three dimensions, Commun. 

Math. Phys., 284 (2008), pp. 459-479. 

[40] S. Viefers, T. H. Hansson, and S. M. Reimann, Bose condensates at high angular mo- 
menta, Phys. Rev. A, 62 (2000), p. 053604. 



STRONGLY CORRELATED PHASES IN RAPIDLY ROTATING BOSE GASES 

CNRS & Laboratoire de Mathematiques UMR 8088, Universite de Cergy-Pontoise, 
Avenue Adolphe Chauvin, 95302 Cergy-Pontoise Cedex, France. 
E-mail address: Mathieu.Lewinamath.cnrs.fr 



Department of Physics, Princeton University, Princeton NJ 08544, USA. 
E-mail address: rseiringSprinceton.edu 



